arXiv: 1509.03310v2 [astro-ph.CO] 21 Feb 2016 


Prepared for submission to JCAP 


Time-delay Cosmography: Increased 
Leverage with Angular Diameter 
Distances 

I. Jee“, E. Komatsu“’^, S. H. Suyu'^’“, D. Huterer^^® 

“Max-Planck-Institut fiir Astrophysik, Karl-Schwarzschild Str. 1, 85741 Garching, Germany 
^Kavli Institute for the Physics and Mathematics of the Universe, Todai Institutes for Ad¬ 
vanced Study, the University of Tokyo, Kashiwa, Japan 277-8583 (Kavli IPMU, WPI) 
'^Institute of Astronomy and Astrophysics, Academia Sinica, P.O. Box 23-141, Taipei 10617, 
Taiwan 

'^Department of Physics, University of Michigan, 450 Ghurch St, Ann Arbor, MI 48109-1040 
^Excellence Gluster Universe, Boltzmannstrasse 2, D-85748 Garching, Germany 

E-mail: ijee@mpa-garching.mpg.de 

Abstract. Strong lensing time-delay systems constrain cosmological parameters via the so- 
called time-delay distance and the angular diameter distance to the lens. In previous studies, 
only the former information was used in forecasting cosmographic constraints. In this paper, 
we show that the cosmological constraints improve significantly when the latter information 
is also included. Specifically, the angular diameter distance plays a crucial role in breaking 
the degeneracy between the curvature of the Universe and the time-varying equation of state 
of dark energy. Using a mock sample of 55 bright quadruple lens systems based on expecta¬ 
tions for ongoing/future imaging surveys, we find that adding the angular diameter distance 
information to the time-delay distance information and the cosmic microwave background 
data of Planck improves the constraint on the constant equation of state by 30%, on the 
time variation in the equation of state by a factor of two, and on the Hubble constant in the 
flat AGDM model by a factor of two. Therefore, previous forecasts for the statistical power 
of time-delay systems were significantly underestimated, i.e., time-delay systems are more 
powerful than previously appreciated. 
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1 Introduction 


The redshift-distance relation constrains cosmological parameters. The use of strong grav¬ 
itational lens systems with time-delay measurements as a cosmological distance probe was 
suggested by Refsdal in 1964 [1], and has been extensively studied to measure the present 
value of the Hubble parameter, Hq (see, e.g. [2]). The physical quantity of interest here is the 
so-called time-delay distance D^f This is a distance-like quantity, given by a combination 
of three angular diameter distances in a strong lens system: 


Dm = (1 + zl) 


Da{EL)Da{ES) 

Da{LS) 


( 1 . 1 ) 


where zl is the redshift of the lens, while Da{EL), Da{ES) and Da{LS) are the angular 
diameter distances from the Earth to the lens, from the Earth to the source, and from the 
lens to the source, respectively. As each of the distance is proportional to the inverse of 
Hq, Dm is proportional to 1/Hq. Ref. [3] has shown that a precise measurement of Hq is 
important for constraining the equation-of-state parameter of dark energy, w. This motivated 
the measurement of time-delay distances in the recent years. However, the constraints weaken 
significantly when w is allowed to vary with time. 

It is, in fact, possible to extract the angular diameter distance to the lens, Da{EL), by 
combining time-delay measurements with the lens stellar velocity dispersion measurements 
[4, 5]. The physics is simple: The time-delay measurement is the mass estimate of the lens 
galaxy, while the velocity dispersion measurement is the potential estimate. By knowing the 
mass and the potential, we can calculate the physical size of the system, thus the system can 
be used as the “ruler”. Comparing the observed separation of lensed images with the length 
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of the ruler, we can estimate Da{EL). In ref. [4], Paraficz and Hjorth demonstrated this 
for a singular isothermal sphere (SIS) lens without mass external to the lens along the line 
of sight. In ref. [5], we have extended the analysis by including the external convergence, 
allowing the mass profile of the lens to follow an arbitrary power-law, and allowing the 
velocity dispersion to be anisotropic. We have found that the main source of uncertainty for 
Da is the unknown (anisotropic) velocity structure, which affects the normalization of the 
potential of the lens galaxy. We have also found that the mass external to the lens along the 
line of sight (external convergence) has no effect on the inferred Da{EL). In this paper, we 
show that the constraints on the cosmological parameters, especially on a time-varying w, 
improve significantly by including Da{EL). 

The structure of the paper is as follows. In section 2, we introduce the mock catalog of 
the strong lens systems we use, and describe the cosmological model along with the fiducial 
cosmological parameters we use. In section 3, we compare the cosmological constraints we 
expect from lenses to those from the other cosmological distance probes, specifically, Type 
la Supernovae (SNe) and Baryon Acoustic Oscillation (BAO). In section 4 we introduce the 
pivot redshift and Wp-Wa parametrization, and show how lensing distances improve the figure 
of merit for the dark energy equation of state. We combine cosmological information from 
different probes to show the constraining power of strong lenses in practice, and conclude in 
section 5. In Appendix A, we show the constraint on the Hubble constant E[q using the two 
lensing distances. In Appendix B, we show the dark energy constraints assuming that the 
universe is spatially flat. 

2 Method 

Each well-modeled time-delay lens system yields two distance (-like) quantities, Da{EL) and 
The uncertainties of Da{EL) and Da* are dominated by the velocity dispersion and the 
external convergence, respectively. In this work we make an assumption that we can measure 
both distances with 5% uncertainty, which requires a few per cent measurement of the spa¬ 
tially resolved velocity dispersion of the lens galaxy, as well as a good understanding of the 
mass distribution along the line-of-sight, that is obtainable by simulations and observations 
of the lens environment [5, 6]. Regarding the lens mass model, the power-law density prohle 
in ref. [5] is widely used due to its ability to fit the imaging data near the image positions. 
The local density profile is well reconstructed with the model if the images are spatially ex¬ 
tended such that information from thousands of intensity pixels can be used. However, ref. 
[7, 8] have pointed out that the information obtained by the lensed images cannot uniquely 
determine the shape of lens mass profile due to the so-called Source-Position Transformation 
(SPT). Specifically, they focused on the degeneracy between composite density profiles and 
a power-law mass profile, and have shown that fixing the shape of the lens mass profile as 
a power law can break the SPT. However, they have also mentioned that these models can 
be distinguished if more information is available: for example, if more than three images 
with time delays are observed, the degeneracy can be broken as the general SPT does not 
conserve the time delay ratios. In ref. [6], the robustness of the measured time-delay distance 
is tested with power law and composite model under the presence of lens kinematics data 
and shown to be nearly independent of the choice of the model. Ref. [9] has shown that the 
so-called Mass-Sheet Transformation (MST), which is a special case of the SPT, scales the 
time delays by the same factor, and thus conserves the time-delay ratio, can bias the mass 
modeling. However, ref. [10] has empirically shown using the Illustris simulation that most 
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of the early-type galaxies with high (cr > 200kms“^) velocity dispersion, which most of the 
lens galaxies are, show nearly power-law behavior. The MST degeneracy can also be broken 
if additional information on the lens galaxy (e.g. velocity dispersion) is obtained. Ref. [11] 
has shown that the existence of substructures in the lens galaxy perturbs the time delay: 
However, the effect of perturbation on the time delays (<1 day) is typically smaller than 
the currently available time-delay measurement uncertainties, and thus both distances are 
mainly determined by the global mass distribution rather than the substructures. Thus we 
claim that this precision measurements on both distances is possible, but only when good 
quality imaging / kinematics data as well as time-delay measurements are available. The 
correlation between the two distances is negligible, because the uncertainties in the velocity 
dispersion and the external convergence are uncorrelated. 

To study the expected cosmological constraints from lenses, we need to specify the 
distribution of lens and source redshifts. We use the catalog of mock lenses in ref. [12] 
to obtain the redshift distribution of time-delay lenses with double and quadruple images 
expected for the Large Synoptic Survey Telescope (LSST) [13, 14]. Although LSST itself is 
expected to find around ten thousand lensed quasars, there are only about 400 systems that 
would have good time delay measurements [15]. To obtain distances from a lens system with 
a reasonable accuracy, a good mass model of the lens galaxy is also required, as both the 
time-delay distance and the angular diameter distance are sensitive to the mass distribution 
of the lens. Ancillary data both in terms of high-resolution imaging and spectroscopy of the 
lens systems are needed for accurate lens mass modeling. Therefore, we select lens systems 
from the mock catalog with the following criteria for acquiring ancillary data with relative 
ease: (1) the quasar image separation is > 1", (2) the third brightest quasar image has an 
Tband magnitude ruj < 21, and (3) the lens galaxy has m* < 22. The criteria on the quasar 
image separation and brightness make it easier to measure the time delays with high precision 
(uncertainty of a few percent). Furthermore, a sufficiently wide quasar image separation is 
required for extracting the Einstein ring of the quasar host galaxy and measuring the lens 
velocity dispersion for mass modeling. The lens galaxy also needs to be of sufficient brightness 
for measuring the lens velocity dispersion. In this work, we focus on quadruply imaged lens 
systems as they provide more information than doubly imaged systems. After applying these 
criteria to the LSST mock lens sample, we obtain 55 quadruple lens systems as the best cases 
of obtaining ancillary data. 

Figure 1 shows the source and the lens redshift distribution of quadruple lenses in 
our mock catalog. For the visualization purpose, only in this figure, the total number of 
detectable lenses is oversampled by a factor of 5 (based on the catalog from ref. [12]) to 
populate the histogram. The expected constraints reported in this paper are derived from 
the actual distribution of the 55 lenses. 

Since these are the bright lens systems, a fraction of these systems will already be 
discovered in the current imaging surveys. In particular, we expect that ~25% of these 
systems will be discovered in the Dark Energy Survey (DES)^ and the Hyper Suprime-Cam 
(HSC; [16])^ Survey. Eurthermore, we expect a few more quadruple lens systems from the 
northern areas of the HSC Survey that are not covered in DES and LSST. Therefore, even 
though we focus here on the LSST sample, our cosmographic predictions are also relevant for 
the upcoming years before the LSST era as new lens systems in the current imaging surveys 
are discovered and monitored. 

^http: //www.darkenergysurvey.org/index.shtml 

^http://www.naoj.org/Projects/HSC/surveyplan.html 
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Figure 1: The distribution of source (zs) and lens (zl) redshifts of quadruply imaged time- 
delay lenses expected for LSST [12], These lens systems are the best ones to obtain ancillary 
data and measure distances with 5% precision. For visualization purposes, the number of 
the lenses is oversampled by a factor of five only in this figure (thus each pixel does not 
necessarily have a multiple-of-5 value). 


We explore constraints on two variations of ACDM model. Both assume a curved 
universe {flk 7^ 0) and an unknown equation of state of dark energy (w ^ — !)• The first 
model assumes that tc is a constant (otcCDM) with the following cosmological parameters: 

9 ^ (orrCDM model). (2T) 


The second model further assumes that w varies in time with w = tco + (1 — a)wa [17] 
(ou;,CDM): 

9 € {^m,^k,WQ,Wa,h} (otc^CDM model). (2-2) 

We choose the fiducial cosmology following Planck 2015 (11^ = 0.308, 11^ = 0, tco = — 1, 
Wa = 0, and h = 0.678) [18]. 

We use the Fisher information matrix (hereafter the Fisher matrix) to calculate the 
constraining power of the cosmological probes. For a data vector d{9) with a set of parameters 
0, the Fisher matrix F is given by 


= E 

a/3 


dda 


Cov 


a/? 


ddp 
89 , ’ 


(2.3) 


where indices a and jd run over the observables, and Cov(d) is the data covariance matrix. 
For lenses, Cov^^ = 6 ai 30 '~‘^, where 6 a^ is the Kronecker delta as we assume no correlation 
between the different lens systems and between the two measured lensing distances 
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and D/\t of each lens system. The uncertainty in each distance Uq is Uq = O.OSdo as we 
assume 5% precision measurements of both distances. The inverse of the Fisher matrix, F“^, 
gives the covariance matrix of the parameters, and the marginalized uncertainty on the z-th 
parameter is calculated as If the normalization of F increases by a factor of n, 

then the normalization of the parameter covariance matrix decreases by a factor of n, thus 
the error bar on each parameter tightens by a factor of ^/n. 


3 Single-probe constraints combined with the Planck distance prior 


In this section we investigate the constraining power of the time-delay lenses expected from 
LSST in section 3.1, and compare the constraints to that of two other cosmological distance 
probes, the BAO data from Baryon Oscillation Spectroscopic Survey (BOSS) data release 
(DR) 11 in section 3.2, and the SNe data from Joint Light-curve Analysis (JLA) in section 
3.3. We examine the future prospects of constraints from BAO and SNe using LSST data 
in section 3.4. We combine each probe with the CMB distance prior calculated from the 
Planck 2015 result [19]. The distance prior is calculated using the shift parameter, i?shiftj 
and the multipole corresponding to the sound horizon at the moment of last scattering, h. 
The definitions of these parameters are 


-f^shift — Da (^*)/c. 


_ Da{z^) 

L — TT , 

rs{z^) 


(3.1) 


where z* = 1089.94 is the redshift of the last scattering surface, and rs(z*) = 144.89 Mpc is 
the size of the sound horizon at redshift The distance prior compresses information in 
the CMB power spectrum relevant for dark energy to two numbers. 


3.1 Time-delay lenses 


We first show how sensitive Da* and D\{EL) are to wq and as a function of zl and zs- 

In figure 2, we show dlnT/dp (where T = {DA,D^t) and p = {wo,'Wa) ). The larger 
the absolute value of dlnT/dp is, the bigger the unmarginalized sensitivity of the distance T 
becomes to the parameter p, with all the other cosmological parameters fixed at the fiducial 
values. Also, equation 2.3 shows that the information is proportional to dT/dp. The higher 
Zl is, the more sensitive D^t becomes to both wq and Wa for a given zs, and vice versa. We 
find that Da is always more sensitive to Wa, when all the other parameters are fixed at the 
fiducial values. 
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2£ for Da, zs for D^t 


ZL 


Figure 2: Logarithmic derivatives of distances T ={DA,DAt) with respect to the cosmo¬ 
logical parameters p =(wo,iOa), as a function of redshift. Left panel: The solid lines show 
T = Da, while the dashed lines show T =D/\t, with various combinations of the lens and 
the source redshift. Each dashed line corresponds to one zl in range [0.2,1.0] in increments 
of 0.2, and shows —din Di\t/dp/0.01 as a function of zs- We only show zs which is higher 
than the highest zl in this range. Right panel: The solid lines show —din DA/dp/0.01 as 
a function of zl. Each dashed line corresponds to one zs in range [1.0,3.0] in increments of 
0.5, and shows —dlnD/\t/dp/0.01 as a function of zl. We only show zl which is lower than 
the lowest zs in this range. Both panels show that Da is always more informative than D/^t 
for constraining Wa (i.e., \dln DA/dwal > \dInD/\t/dwa\), and Da is often more informative 
than D^t on wq. 




Figure 3: The 1-cr uncertainty in w, denoted as cr(w), and that in wg, denoted as (T(t(;o), 
from time-delay lenses as a function of the number of lenses for the left and the right panel, 
respectively. The black dash-dot line is the lens-only data, while the other lines use the 
Planck distance priors combined with Da (blue dashed), D^t (red dotted), or both (green) 
from lenses. (Left) otcCDM model. (Right) orCzCDM model marginalized over Wa, as well 
as all the other parameters. 

In figure 3, we show the expected l-cr uncertainties in w for the otcCDM model and wg 
for the orc^jCDM model (combined with the Planck distance priors), with all the other pa- 
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rameters marginalized over. As the Fisher matrix is proportional to nqu^ds) the marginalized 
uncertainty in Da{EL) + D^t scales as oc l/^Uquads- For otcCDM (left panel of figure 3), 
the Planck distance priors combined with either Da (blue dashed line) or (red dotted) 
from lenses improve the constraint on w significantly compared to the lens-only case (black 
dash-dot). Combining all improves the constraint further by 30%, in comparison to that of 
the Da* + Planck. For otc^CDM (right panel of figure 3), we find that the lens-only {Da + 
Dm) breaks the degeneracy between wq and the other parameters more efficiently than either 
combination of Planck -I- Da or Planck -I- Dm-, yielding a tighter constraint. Combining all 
improves the constraint further by a factor of 2 in comparison to Dm + Planck. 




Figure 4: (Left) The 68 per cent CL in the WQ-Wa plane constrained from 55 lenses for the 
otc^jCDM model. The blue line is the marginalized constraint from Da, the red line is the 
marginalized constraint from Dm, and the green line is the marginalized constraint from the 
combination of the Da and Dm- The black line is the unmarginalized constraint from Planck 
only. The horizontal and vertical dashed lines correspond to the fiducial values of rco = — 1 
and Wa = 0, respectively. (Right) Zoom in of the left panel, with degeneracy directions 
plotted as arrows. The arrows show the response of wo and Wa to the shift AOfc = 0.02. The 
blue arrow is the response for Da{EL), the red arrow is for Dm, and the green arrow is for 
the combination of two. The larger the arrow is relative to the contour, the more sensitive 
the distance is to The alignment between the arrows and the major axis of the contours 
indicate that the constraints from Da (the red arrow and the red contour) and Dm (the blue 
arrow and the blue contour) are individually dominated by 12^. However, the degeneracy 
between and w is shown to be broken when Da and Dm are combined. For visualization 
purposes, the sizes of the arrows are inflated by a factor of 100. 

The left panel of figure 4 is useful for understanding these results. The marginalized 
uncertainties in Wa from either Da or Dm individually are similar. The Planck distance prior 
(the black line) provides a degenerate combination of wq and Wa, thus cannot be marginalized. 
However, it is nicely orthogonal to the ones from Da (the blue contour) and Dm (the red 
contour). Thus, including the combination of Da and Planck with the previous lensing 
constraints from Dm reduces the uncertainty in Wa significantly. We also note that the 
constraint on Wa from Da + Dm (the green contour) is significantly tighter than the naive 
addition of the blue and the red contour, which indicates that this combination of distances 
effectively breaks the degeneracy between the equation of state of dark energy and the other 
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parameters over which we marginalize. 

Next, we study the degeneracy structure of parameters constrained from lensing dis¬ 
tances, by shifting a parameter and calculating the response of the other parameters to the 
shift in order for the likelihood to be maximized. Specifically, we describe the degeneracy 
between the curvature density Qk and the equation of state parameters wq and Wa- When 
is shifted by an amount to maximize the likelihood at the new fiducial value 
all the other parameters have to be shifted accordingly. The general expression for the shift 
in an arbitrary parameter A0j due to a shift in a fixed, single parameter AOf^ that maximizes 
the likelihood can be calculated as 

= (3.2) 

where is the inverse of the Fisher matrix. In figure 4, the right panel shows the projection 
of the shift vectors to WQ-Wa plane as arrows, along with the marginalized constraint contour 
at the fiducial parameter to display it quantitatively. For each of Da and D^f, the degeneracy 
directions are parallel to the major axes of the contours, which indicates that the degeneracies 
with rifc dominate the dependences of wq and Wa to other parameters. However, when the 
two distances are combined, the curvature degeneracy breaks and the alignment between 
the error contour and the shift disappears (the green arrow and the green contour are not 
aligned). The relative size of the arrow to the contour shows the sensitivity of the probe 
to Hfc: the bigger the vector is with respect to the contour, the more sensitive the probe is 
to the change in By comparing the relative size of the red and the blue arrows to the 
red and the blue contours, we show that Da and Da* are comparably sensitive to AH^, but 
the combination of two increases the sensitivity significantly (the green arrow and the green 
contour). 






Wo 


Figure 5: The marginalized 68 per cent CL constraints from 55 lenses in the (left) Vtm-w 
plane for the orcCDM model, and (right) WQ-Wa plane for the orc^CDM model. The black 
dash-dot lines show the lens-only constraints from Da + D^t, the blue dashed lines the 
constraints from Da + Planck, the red dotted lines the constraints from Da* -|-Planck, and 
the green solid lines the combination of the two distances -|- Planck. 


The left panel of figure 5 shows the joint constraints on and w for orcCDM. We find 
that the Planck distance prior plays an important role in constraining while using both 
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Da and Da* combined with Planck distance prior improves the constraint on w by about 
30% compared to the case of Di^t combined with Planck distance prior. The right panel of 
figure 5 shows the same for wq and Wa for ouJzCDM. We also show the lensing distances + 
Planck constraint on Hq in Appendix A. 

Next, we compare these constraints with those from Planck + BAO and Planck + SNe. 
We calculate the constraints from BAO and SNe using currently available data (BOSS DRll 
for BAO, JLA sample for SNe). 

3.2 BAO 

BOSS DRll provides the volume-averaged distance, Dy = {cz{l + z)‘^D\/H^ ' , at two ef¬ 
fective redshifts (0.32, 0.57) obtained from the BAO peak position in the spherically averaged 
two-point functions. The lower redshift is the LOWZ sample, and the higher redshift is the 
CMASS sample. Also, by separately measuring the two point functions along the line of sight 
and the direction perpendicular to it, the DRll CMASS sample separately constrains the 
angular diameter distance Da and the Hubble parameter 77 at z = 0.57 [20, 21]. To account 
for the correlation among Dy, Da and H, we use the full likelihood of the CMASS sample 
for the analysis provided by the BOSS collaboration [20]. 

In our analysis we assume that the sound horizon scale at the baryon drag epoch, rg^drag, 
is fixed as Vg^drag = 149.28 Mpc [20]. We then combine the cosmological constraints from Dy 
at z = 0.32, and Da and H at z = 0.57 with the Planck distance prior. We calculate the 
Fisher matrix by taking the derivatives of the log likelihood at the fiducial cosmology. The 
results are shown in figure 6. The precision of the BAO data yields the narrowest contours 
on the Vlm-w (for ou^CDM) and WQ-Wa (for ouj^CDM) planes. However, due to the limited 
number of redshifts (z = 0.35 and z = 0.57), the degeneracy is not broken efficiently; thus, 
the expected Planck -|- lens from 55 lenses can improve the constraints significantly, even 
though the precision of lensing data per redshift is not as precise as BAO. 

3.3 SNe 

We now study the constraints from Planck -|- SNe. We use the JLA data [22] to calculate the 
constraints from SNe. JLA uses Supernovae Legacy Survey (SNLS), Sloan Digital Sky Suvey- 
H (SDSS-H) Supernova survey and a few low-redshift samples. The redshift of subsamples 
are: the low-redshift sample (z < 0.1), SDSS-H (0.05 < z < 0.4), and SNLS (0.2 < z < 1). 
There are 740 spectroscopically confirmed type la SNe in JLA. SDSS-H is used for anchoring 
the distances, and also an empirical relation between the host galaxies and the supernovae 
brightness is used as an extra calibration for the absolute magnitude of the SNe. For the 
calibration, there are 4 additional nuisance parameters that are taken into account in JLA: 
a, which scales the stretch of the light curve in time-domain; /3, which scales the color at the 
peak of the light curve; M, which is the absolute B band magnitude of the SNe at the peak 
of the light curve; and Am, which characterizes the peak absolute magnitude change due the 
stellar mass of the host galaxy. 

We use Montepython [23] to sample the JLA likelihood. Specifically, we run Markov 
Chain Monte Carlo to sample the likelihood surface, and compute the covariance matrix in 
the cosmological parameters. We then use its inverse as the JLA Fisher matrix. The results 
are shown in figure 6. While the absolute distances, such as those from BAO and lenses, are 
effective at measuring when combined with CMB [24], the relative distances from SNe are 
not. Thus, when is set free, the constraints on for orcCDM from 55 lenses combined 
with Planck are significantly better than those from 740 SNe combined with Planck. 
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Figure 6: The marginalized 68 per cent CL constraints from strong lenses, SNe, and BAO, 
each combined with the Planck distance prior, in the (left) plane for the orcCDM 

model, and (right) WQ-Wa for the orCzCDM model. The green lines show the constraints from 
Planck + lens, the red lines Planck + JLA, and the black lines Planck + BOSS. 

3.4 Comparison to future BAO and SNe predictions 

With several billion galaxies expected to be detected with LSST, BAO will allow measure¬ 
ments of distances with ~2% precision in the redshift range 1 < z < 3 [25]. Combined with 
Planck, BAO will constrain wq with uncertainty ~0.4 and Wa with ~1 [26]. Also, 500,000 SNe 
are expected to be detected in 10 years of LSST operation in the redshift range 0.1 < z < 1.2. 
With a subsample of 50,000 SNe only, the data will constrain wq with uncertainty ~0.05, and 
Wa to order unity, assuming a flat universe [13]; in combination with Planck, the full sample 
of SNe constraints will be ~0.25 for wq and ~1.2 for Wa for the orc^CDM model [26]. We note 
that a modest sample of 55 lenses combined with the Planck distance prior constrains wq 
and Wa to ~0.4 and ~1.2, respectively (see, e.g.. Figure 6), which is comparable in precision 
to those expected from future BAO or SNe samples in the LSST era, when each is combined 
with Planck. Therefore, strong lenses provide an independent and competitive probe of dark 
energy. Needless to say, leasing, SNe, and BAO are affected by different systematic errors, 
and thus cross-checking the results using these three low-redshift probes of the expansion of 
the universe will be powerful. 


4 Pivot redshift 

The equation of state of dark energy, w{z), can be re-written as 

W = Wq + {1 - a)Wa = Wp + {ttp - a)Wa, (4.1) 

where Wp = wq + {1 — ap)wa [27]. In this parameterization, the pivot redshift Zp = -^— 1 
is defined as the redshift where the uncertainty in w is minimized. The uncertainty in Wp 
shows how well a probe can measure the equation of state w, as Wp is orthogonal to Wa by 
construction, and thus is not coupled to the time variation of w [27, 28]. Zp shows at which 
redshift the main constraint on w is coming from: this pivot redshift varies depending on the 
probes, their redshift distributions and the measurement uncertainties, and can be negative. 


- 10 - 




















Cosmological model 

Probe 

Zp 

(T{Wp) 

(y{Wa) 

FoM 

Flat, fixed flm 

CMB+BAO+SNe 

0.288 

0.0375 

0.364 

73.3 

Lens+CMB+BAO+SNe 

0.487 

0.0296 

0.166 

204 

Flat, marginalized over 0^ 

CMB+BAO+SNe 

0.358 

0.0487 

0.413 

49.7 

Lens+CMB+BAO+SNe 

0.386 

0.0310 

0.363 

88.9 

Marginalized over Om and 

CMB+BAO+SNe 

0.215 

0.0625 

1.03 

15.5 

(ow^CBM) 

Lens+CMB+BAO+SNe 

0.245 

0.0479 

0.621 

33.6 


Table 1: Pivot redshift, Zp, uncertainties in Wp and Wa, and the Figure of Merit (FoM) for 
Planck + BOSS + JLA and Lens + Planck + BOSS + JLA. We test three cases: (1) 0^ = 0 
and = 0.309, (2) = 0 and marginalized over and (3) marginalized over both 

and rim- As there is no correlation between Wp and Wa, FoM = 1 /{a{wp)a{wa))■ 



Wp 

Figure 7: Same as figure 4, but in Wp-Wa plane. Zp is calculated separately for each probe: 
Zp{Da) = -0.226, Zp{DAt) = 0.0251, and Zp{DA + Dai) = 0.0759. 


Figure 7 shows the constraints on {wp,Wa) using the lens distances alone in Wp-Wa 
plane. The constraint on Wp from Da alone is much weaker than that from Dai, while 
the constraints on Wa from both distances are comparable. However, by combining the two 
distance measures, the constraint on Wp improves by a factor of 2, and that on Wa improves 
by a factor of 24, due to complementary degeneracy directions as shown in figure 4. 

The pivot redshift, Zp, and the uncertainties in Wp and Wa for two combinations of 
probes (Planck + BOSS + JLA, and 55 lenses (denoted by Lens) + Planck + BOSS + 
JLA) for three different cosmological cases ((1) = 0 and fJm = 0.309, (2) 0^. = 0 and 

marginalized over Om, and (3) marginalized over both and flm) are summarized in table 
1. In comparison to Planck + BOSS + JLA, adding 55 lenses tightens the FoM by a factor 
of 2.78, 1.62 and 2.32 for the three models, respectively. In particular, when is allowed to 
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vary, the uncertainty in Wa is reduced by almost a factor of 2 by including the lenses, which 
shows that the combination of Da and -Da* is powerful in breaking the degeneracy. 

This is consistent with our argument in section 3.1. Also Zp becomes higher as we include 
the lens distances in every case, which is typically beneficial when combining these probes of 
geometry with probes of the growth of cosmic structure. 


5 Conclusion 


There is more valuable cosmological information in the strongly lensed systems than mea¬ 
surements of the Hubble constant from time delays. In this paper, we have demonstrated 
that the addition of the angular diameter distance measurements to the quantity {D/^t) that 
captures the cosmological information from time delays in the same sample of lenses provides 
crucial help in breaking cosmological parameter degeneracies. This improvement is most sig¬ 
nificant in some of the most interesting parameterizations that are currently being studied, 
such as when curvature of the universe and the time-variation in the equation of state of 
dark energy are allowed to be nonzero. 

We have calculated the lensing constraints based on the predictions for the LSST survey, 
adopting a catalog of 55 quadruply imaged lenses (out of a much larger total number) that 
should have sufficiently good information that all observable quantities of interest in the 
lenses can be accurately measured. We have combined the forecasted lensing information 
from both the angular diameter distance and the time-delay distance. We then compared 
this lensing constraints with that from the BOSS DRll and from the JLA type la supernova 
sample - each combined with the Planck 2015 distance prior. 

We find that the combined lensing information significantly helps constrain the cosmo¬ 
logical parameters, particularly when curvature is allowed to vary and when the equation 
of state of dark energy is allowed to be time-dependent. For example, lensing information 
would improve the current BAO-|-CMB-|-SN constraints on Wa by a factor of two, and those 
on the overall figure of merit of dark energy by about a factor of two relative to the case with 
no lensing (see Table 1 and figure 8). Key to this significant improvement is lensing’s ability 
to break the degeneracy between curvature and the equation of state parameters; see figure 
4. 


We are therefore very optimistic about the prospects of a select, accurately observed 
subsample of strong gravitational lenses to improve our constraints on dark energy. Fortu¬ 
nately, the lensing samples are a guaranteed product of the current and upcoming wide-held, 
deep surveys such as HSC, DES and LSST. 
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Figure 8: The marginalized 68 per cent CL from strong lenses, combined with Planck, BAO 
and SNe in the (left) VLm-w plane for the orcCDM model, (right) WQ-Wa for the otCzCDM 
model. 
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A Lensing constraints on Hq 

A.l Hq in ozuCDM and ou^xCDM models 

We show the constraints on Hq from Da and D^t combined with the Planck distance prior. 



Figure 9: The l-cr uncertainty in h = Hq/100 km/s/Mpc from time-delay lenses as a 
function of the number of quadruply imaged lenses for the (left) orcCDM model, and (right) 
otc^jCDM model. 

Figure 9 shows the expected 1-a uncertainties in h from strong lenses combined with 
the Planck distance priors. As D^^t is mostly sensitive to Ifo, the constraining power of 
D^t + Planck (red dotted line) is more powerful than that of Da + Planck (blue dashed 
line). When w is fixed as a constant (orcCDM model, left panel), Ha* + Planck are more 
powerful than D/\t + Da (black dot-dashed line). When w is allowed to vary (otc^CDM 
model, right panel), however, Da* + Da is more powerful than Da* + Planck. This is due 
to the degeneracies between Ho, flk and w from the linear CMB constraints alone [29], which 
cannot be broken by I?At- However, ref. [3] has shown that the main degeneracy from CMB 
constraints is between w and Hq, and as shown in section 3.1, the combination of lensing 
distances is powerful in breaking the degeneracy between and w. Thus, the combination 
of Planck and the lensing distances shows 30% improvement in constraining h. 
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A.2 Hq in flat ACDM model 

We show the constraints on Hq for the ACDM model in figure 10. 



Figure 10: Same as figure 9, but for the flat ACDM model. Here the blue dashed, red dotted 
and black dash-dot lines are from the leasing distances alone, not combined with Planck. We 
show the constraints from Planck as the horizontal magenta solid line and Planck -|- leasing 
distances as the green solid line. Planck-precision constraint in h is achievable with 10 lenses 
when we use both Da and Haj, while we need 25 lenses to achieve the same constraint from 
D^t alone. 


Assuming that 5% precision measurements in individual distance (both Da and Haj) 
are achievable from lens systems, 10 lenses are enough to measure the Hubble constant to 
the same precision as Planck. The number of required lenses to achieve the same precision 
increases to 25 if constraints are from D^^t only. 


B Constraints assnming the flat universe 

In section 3.1, we have shown that the leasing distances are powerful probes for the curvature 
of the universe. Specifically, Da and D^^t respond to curvature differently in WQ-Wa plane, 
thus the combination of two gives a strong constraint on We repeat the same analysis 
for the flat universe model (Dfc = 0). The model parameters are summarized as 

(flat tcCDM model), (B-1) 

and 

0 G {Dm, tco, rca, (flat rCzCDM model). (B.2) 

The constraining contours for these models are shown in figures 11 and 12. 
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Figure 11: Same as figure 5, but for the flat (left) wCDM and (right) tu^CDM model. 




Figure 12: Same as figure 6, but for the flat (left) ruCDM and (right) u^xCDM model. 


Figure 11 shows that under the flatness assumption, the constraints from Dj\t + Planck 
are already as tight as those from Da + + Planck, i.e., the constraining power from 

Da in flat universe is minor. Figure 12 shows that the 55 lenses combined with Planck still 
constrain the equation of state better as compared to Planck + JLA and Planck + BOSS 
for the flat rcCDM model (left panel), and comparably well as Planck + BOSS for the flat 
tc^CDM model (right panel). 
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